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ABSTRACT 

We report on the discovery of a surprising observed correlation between the slope of the low-mass stellar 
global mass function (GMF) of globular clusters (GCs) and their central concentration parameter c = log(r t /r c ), 
i.e. the logarithmic ratio of tidal and core radii. This result is based on the analysis of a sample of twenty 
Galactic GCs with solid GMF measurements from deep HST or VLT data. All the high-concentration clus- 
ters in the sample have a steep GMF, most likely reflecting their initial mass function (IMF). Conversely, 
low-concentration clusters tend to have a flatter GMF implying that they have lost many stars via evaporation 
or tidal stripping. No GCs are found with a flat GMF and high central concentration. This finding appears 
counter-intuitive, since the same two-body relaxation mechanism that causes stars to evaporate and the cluster 
to eventually dissolve should also lead to higher central density and possibly core-collapse. Therefore, more 
concentrated clusters should have lost proportionately more stars and have a shallower GMF than low concen- 
tration clusters, contrary to what is observed. It is possible that severely depleted GCs have also undergone 
core collapse and have already recovered a normal radial density profile. It is, however, more likely that GCs 
with a flat GMF have a much denser and smaller core than suggested by their surface brightness profile and 
may well be undergoing collapse at present. In either case, we may have so far seriously underestimated the 
number of post core-collapse clusters and many may be lurking in the Milky Way. 
Subject headings: Stars: luminosity function, mass function — Galaxy: globular clusters: general 



1. INTRODUCTION 

In recent years an increasing number of globular clus- 
ters (GCs) have been found to be seriously depleted 
in low-mass stars (< 0.5 M ffl ) when compared with 
high-conce ntratio n c lusters dParesce & De Marchil I200C 
herea fter iPDMt iDe Marchi. Paresc e & Portegies Zwart 
2005). The first heavily depleted cluster t o be discovered 
was NGC 67 12 dDe Marchi etal.1 1 19991: lAndreuzzi et al.l 
12001. followed by Pal 5 (iKoch et all 120041). NGC 6218 
dDe Marchi, Pulone & Parescd 120061 h ereafter |PMPP), 



2007) and NGC6838 
In all cases, the 



NGC 2298 dDe Marchi & Pulond . 
dPulone. De Marchi & Parescd I2007F 
analysis of the radial variation of the directly observed 
mass function (MF) confirms that the relative depletion of 
low-mass stars is not due to the local effect of mass segrega- 
tion but is a structural property of the global mass function 
(GMF), i.e. the MF of the entire stellar population of the GC 
obtained by model fitting. Hereafter, the GMF is defined as 
a power-law function, where the number of stars N per unit 
mass m follows a relationship of the type dN /dm oc m a for 
0.3<m<0.8M Q . 

For some of the clusters (NGC 6712, Pal 5, NGC 6838), 
this finding is at least qualitatively consistent with the 
predicted effect of tidal stripping caused by the Galaxy 
as these objects have typically shorter disruption times 
( J dis) than the average GC according to the models 
of iGnedin & Ostrikerl d 19971). iPinescu. Girard & van Altenal 
d!999l) and iBaumgardt & Makinol (120031) . Quantitatively, 
however, the agreement is poor for all of them, since 
the expected tight correlatio n between 7di s and GMF slope 
(Baumgardt & Makino 2003), whereby clusters with a higher 
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probability of disruption sho uld alw a ys have a shallower 
GMF slope, is not observed (D MPPl foe Marchi & Pulonel 
120071) . Rather than an error in the models, this is possibly 
the result of the large uncertainties affecting their input pa- 
rameters, especially the clusters' orbits and exact shape of 
the Galactic potential. (For example, the redetermination of 
the orbit of NGC 62 18 with respect to the Hipparcos refer- 
ence system has led to a value o f Zdi s in better agreement with 
the cluster's flat GMF (DMPP) than that based on the previ- 
ous, less accurate orbit.) This situation is unlikely to improve 
much until the advent of space astrometry missions like Gaia. 

The absence of a clear correlation with the effects of tidal 
stripping, on the other hand, could also be due, at least in 
part, to our imperfect understanding of the relation between 
the cluster's GMF and its fundamental structural parameters 
and their evolution in time. To explore this possibility, we 
have used the available data to search for possible signs of a 
more complex situation. We present in this Letter the results 
of this study which indeed imply that there is more at work 
here than was thought up to now. 

Specifically, we find that there is a systematic trend be- 
tween the GMF slope a and the central concentration pa- 
rameter c (defined as log(r t /r c ) where r t and r c are the 
tidal radius and core radius, respectively), in the sense that 
all five clusters above with a severely depleted GMF have 
an intermediate or low value of c. On the co ntrary, the 
twelve halo clusters in the sample studied by PDM| and 
IDe Marchi. Paresce & Portegies Zwartl d2005l) have typically 
high central concentration (<c> ~ 1.9) and correspondingly 
much steeper GMF slopes, typically a = -1.4 in the range 
0.3 — 0.8 M Q . In the following sections, we explore in more 
detail the origin and nature of this trend and its possible ex- 
planations. 

2. THE SAMPLE 

The sample used in this investigation includes all GCs for 
which reliable luminosity functions (LFs) and MFs exist from 
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TABLE 1 

Mass function index and central concentration 



Object Reference a c My 

-9.42 
-6.74 
-6.30 
-5.17 
-10.29 
-8.93 
-7.20 
-7.32 
-7.48 
-8.20 
-6.48 
-6.63 
-7.23 
-8.50 
-7.50 
-7.73 
-7.55 
-5.60 
-9.17 
-7.43 



NOTE. — For clusters marked with an asterisk following the cluster 
name the value of a is that of the GMF. For all other objects, a is the 
index of the MF measured near the half- mass radius. Bibliographical ref- 
erences in Column (2 ) ar e as follows: (a) Paresce & De Marchi 2000; (b) 
IBellazzini et"aT]|2002l; (c) iDeM archi & Pulone 2007; (d) Ko ch et alj|2004t 
(e) Pulone, De Marchi & Paresce 1999; (f) De Marchi, Pulone & Paresce 
120061; (s) IPulone et alj I2003t (h) IDe Marchi. Paresce & Pulone! I2000t 
(T) lAlb row. De Marchi & Sahu 200 2); (i) lAndreuzzi et al.1 120011 ; (k) 
IPulone. De Marchi & ParesceH2007t (p fPasquafi et aU2004l 

deep HST or VLT photometry. This includes the twelve halo 
clusters studied by PDM, the two bulge cluste rs NGC 6352 
and NGC 6496 studied by IPulone et al.1 (120031) and the five 
GCs, mentioned above, that have recently been shown to 
have a depleted GMF at the low-mass end. Finally, we have 
ad ded to the sampl e NGC 288, whose LF has been studied 
by IBellazzini et al.1 (120021) with the HST and also reveals a 
paucity of low-mass stars. The complete sample is listed in 
TableQ] where column (1) gives the cluster name, column (2) 
the bibliographic reference, column (3) the MF index a in the 
mass range 0.3- 0.8 M Q , column (4) the value of the central 
concentration parameter c and column (5) the total integrated 
magnitude My, both from Harris (1996). 

As for the value of a, it has been derived as follows. For 
those clusters for which a GMF index exists from multi-mass 
models (indicated by an asterisk following the cluster name), 
that value is used. For all other clusters, a is the index of 
the power-law MF that, once folded through the derivative of 
the mass-luminosity (M-L) relationship appropriate for the 
cluster's metallicity, best fits the LF measured near the half- 
mass radius (rhm), since there the MF is expected to cl osely 
approach the GMF dDe Marchi. Paresce & Pulondl2000l) . For 
clusters in the sample studied bv lPDML we adopted the same 
M-L used in that w ork, while for the re maining objects the 
M-L relationship of Bar affe et alj ( 1 19971) for th e appropriate 
metallicity was used. As pointed out by iPDMl since the LF 
of NGC 6341 and NGC 7099 were measured farther out in the 
cluster, namely at about four times the half-mass radius, in 
principle a small positive correction to the measured index a 
should be applied to account for the steepening of the MF 
with increasing radial distance caused by mass segregation. 
This correction amounts to less than 0.2 dex and is included 
in the values listed in TableQ] The correction being relatively 
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FIG. 1. — Observed trend between MF index a and the central concentra- 
tion parameter c. Clusters are indicated by their NGC (or Pal) index number. 
Objects for which a GMF index is available are marked with a circled cross. 
The large crosses indicate the average and 1 a distribution of a and c for clus- 
ters with c values above and below 1 .4. The dashed line is an eye-ball fit to 
the distribution. 

small, however, none of our conclusions would change if we 
were to ignore it. 

In FigureQ] we show the run of a as a function of c from 
TableQ] It is immediately obvious that there are no high- 
concentration clusters with a shallow GMF. It also appears 
that a relatively low concentration is a necessary but not suf- 
ficient condition for a depleted GMF. The median value of 
c (1.4) splits the cluster population roughly in two groups, 
one with lower and one with higher concentration. The mean 
GMF index of the first group is a = -0.3 ± 0.7, while the 
second has a much tighter distribution with a = -1.4 ±0.3. 
The average values and the associated ± 1 a uncertainties are 
shown as thick crosses in FigureQ] Due to the large disper- 
sion of a at low concentration, it is not possible to derive 
a precise correlation between c and a over the whole range 
spanned by these parameters. On the other hand, the GMF 
becomes undoubtedly less steep as c decreases. The relation- 
ship a + 2.5 = 2.3/c is a simple yet satisfactory eye-ball fit to 
the distribution (dashed line in FigureQ]i. 

3. DISCUSSION 

The result shown in Figure 1 is very surprising, since there 
is a clear absence of centrally concentrated clusters that are 
depleted in low-mass stars as one would expect of GCs un- 
dergoing tidal disruption, i.e. objects that should occupy the 
upper right quadrant of FigureQ] This is contrary to theoret- 
ical expectations, since the same relaxation mechanism that 
drives a cluster towards higher central density and, eventu- 
ally, core collapse is als o responsible f or the dissolution of the 
cluster, via evaporation (Spitz erTl987l) . Therefore, as a cluster 
evolves dynamically in the course of its lifetime, one would 
expect that its GMF should become shallower at the low- 
mass end while the central concentration parameter should 
increase. Put differently, severe mass loss should be a suf- 
ficient, yet not necessary, condition for core collapse. In this 
simple scheme, one would for instance expect that NGC 2298 
and NGC 6218, with a ~ 0.3-0.5 should have a much denser 
core than NGC 6397 with a ~ -1.6. Clearly, this picture is 
not at all consistent with the results shown in FigureQ] 

That cluster concentration somehow creates variations in 
the IMF is not a plausible explanation for the observed trend 
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Why haven't loose globular clusters collapsed yet? 
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between c and a, since the star formation process could not 
possibly know about the future structural properties of the 
forming cluster. Thus, unless the high c value of the densest 
clusters is primordial and a presently unknown mechanism 
exists whereby low-mass star formation is hampered in low 
density environments, the depletion at the low-mass end of 
the GMF must result from mass loss via dynamical evolution. 

It is possible that the clusters that have lost a significant 
fraction of their original mass have indeed undergone core 
collapse and have already recovered a normal radial den- 
sity and surface brightness profile. A large body of theo- 
retical studies exists on the evolution of GCs after core col- 
lapse dHenonll 1975!; ISugimoto & Bettwieserill983t iGoodmanl 



1993; Makinol ll996i r according to which they should un- 
dergo a homologous re-expansion, within a thermal times cale, 
triggered by the energy released by hardening binaries dHuj 
1985). This process, however, is not necessarily stable 
and can lea d to re-collapse and subsequent gravothermal 
oscillations dSu gimoto & Bettwieser 1983; Goodman 119871; 
ICohn, Hut & Wisdll989tlMakinolll996l) when the number of 
stars in the clust er exceeds a threshold that d epends on the 
mass spectrum. iMurphv. Cohn & Hutl fl990) show that for 
a Salpeter IMF the re-expansion of the core is stable up to 
iV s ~3 x 10 5 , a value comparable with the estimated number 
of objects present in the most depleted clusters in our sam- 
ple. It, therefore, appears plausible that these clusters may 
have re-expanded after core collapse. On the other hand, 
Murp h~Cohn & Hull (fl99 0) predict the observed core radius 
to expand over time as r c oc f 6 and it would take too long for 
the core to reach a size comparable to that of the pre-collapse 
phase, unless the core shrunk only marginally during collapse 
or current models under-estimate the observed size of r c dur- 
ing the collapse phase of a realistic multi-mass cluster with a 
finite number of stars. 

The evaporation timescale increases with increasing total 
cluster mass Mt, so one could possibly understand the trend 
shown in FigureQ] if clusters of lower central concentration 
were also those with lower mass. Unfortunately, this is not 
easy to investigate as less than half the clusters in our sample 
have a relatively solid Mj, based on multi-mass model fitting, 
while for the rest Mt comes from the total luminosity, i.e. the 
integrated magnitude My, by assuming a constant mass-to- 
light ratio. The value of My is shown for each cluster in Ta- 
bleland no correlation is found between it and c. Similarly, 
no correlation exists between Mj and c for clusters with a re- 
liable value of M T . For example, NGC 2298, NGC 6218 and 
NGC6656 have masses in the ratio 1 : 2.4 : 5.3, but share the 
same concentration c ~ 1.3. Conversely, NGC 6397 is about 
half as massive as NGC 6656, but its c is twice as large. 

It may be argued that besides evaporation other tidal mech- 
anisms are responsible for the depleted GMFs that we find. 
In fact, the compressive heating that GCs undergo when 
they cross the Galactic plane (disk shocking) or venture 
close to the Galactic centre (bulge shocking) can have a 
much stron ger effect than evaporation, depending on the clus- 
ter's orbit dAguilar. Hut & Ostrike"ril 19881: [Gnedin & Ostrikeri 
1 19971; iDinescu. G irard & van AltenaT ll999l) . However, while 
bulge and disk shocks can cause significant mass loss 
dSpitzerl 119871; iHeggie & Hull 120031) . they will preferen- 
tially remove low-mass stars only if these have pre- 
viously been push ed toward the cluster's p eriphery by 
mass segregation ( Ves perini & Heggiel 119971) . Even if 
for some GCs tidal shocks have been as important as 
(or more important than) internal two-body relaxation in 



determining the mass loss rate, the observed trend is 
still puzzling since in any case tidal shocks should ac- 
celerate the evolution of a cluster t oward higher cen- 
tral density and core collapse (see ISpitzer & Chevalier! 
119731: IChernoff. Kochanek & Shapirol 119861: ISpitzerl 119871) . 
Only initially very loose GCs (c < 0.5) are expected to 
quickly dissolve, via mass lo ss due to stellar evolution , 
before reaching core collapse (IChernoff & Weinberg! |l99"ot 
iFukushige & Heggielll995l) . 

However, FigureQ] shows that also clusters of interme- 
diate concentration may undergo severe mass loss without 
necessarily showing signs of core collapse for as long as 
a Hubble time. This behaviour is in principle consistent 
with the results of Fokker-Planck and N-body calculations 
of realistic clusters, includin g the effect of stellar evolu- 
tion and two-body relaxation. [Chernoff & Weinberg! (1 1990!) 



and lTakahashT& Portegi es Zwartl ( 20001) have investigated 
the temporal evolution of c and a for various initial condi- 
tions (total mass, concentration, relaxation time and IMF in- 
dex). Their models suggest that, unless the IMF is very steep 
(a ~ -3.5), energy equipartition and mass segregation will 
initially drive the cluster towards lower values of the cen- 
tral density, mainly because the cluster shrinks (and therefore 
stars are lost) more quickly than the core can contract. Even- 
tually the cluster undergoes core collapse, but how long this 
takes and whether the increase in the central concentration c 
can measurably affect the surface brightness profile depends 
on the initial conditions and in particular on the IMF slope. If 
the latter is very shallow (a ~ -1 .5), stellar evolution may re- 
move enough mass from the cluster so that its core is reduced 
to just a few stars. Even in the deepest collapse phase, no 
central cusp would be vi sible in the surface brig ht ness p rofile. 

The predictions of IChernoff & Weinberg! (1 19901) and 
Takah ashi & Portegies Zwartl (l2000h are difficult to compare 
directly to the actual data, since the measured value of c, 
based on the surface brightness profile dominated by red-giant 
stars, is not a good tracer of the true cluster density. Neverthe- 
less, the overall picture seems compatible with FigureQ] and 
their findings allow us to put forth the following hypotheses 
to explain the presence of clusters with a shallow GMF and 
low c and the absence of objects with a shallow GMF and 
high c. 

The dashed line in Figure 1 approximately traces the evo- 
lutionary path of GCs from their birth towards two opposite 
directions of increasing or decreasing concentration. Clus- 
ters born with sufficiently high concentration (c > 1 .5) evolve 
towards core collapse. Mass loss can be important via stel- 
lar evolution in the first ~ 1 Gyr, and to a lesser extent via 
evaporation or tidal stripping throughout the life of the clus- 
ter, but the GMF at any time does not depart significantly 
from the IMF. Clusters with c < 1 .5 at birth also evolve to- 
wards core collapse, but mass loss via stellar evolution and, 
most importantly, via relaxation and tidal stripping proceeds 
faster, particularly if the orbit has a short perigalactic distance 
or frequent disk crossings. Therefore, as the tidal boundary 
shrinks and the cluster loses preferentially low-mass stars, the 
GMF progressively flattens. This speeds up energy equiparti- 
tion, but c still decreases, since the tidal radius shrinks more 
quickly than the luminous core radius (although the central 
density, particularly that of heavy remnants, is increasing). 
These clusters can eventually undergo core collapse, but this 
might only affect a few stars in the core, thereby making it 
observationally hard to detect. 

An alternative possibility is that there were IMF variations. 
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Some clusters with shallow IMF have undergone severe stel- 
lar mass-loss and have therefore expanded considerably. This 
has led to a lower c and a shallower GMF (larger a) because 
these systems were more prone to tidal stripping. Most of 
these clusters have already disrupted but some survive for a 
long time in a state of low c and large a before collapse oc- 
curs. Clusters with a steeper IMF, on the other hand, have 
proceeded normally to core collapse. 

The problem with this alternative possibility is the absence 
of GCs with a shallow GMF slope and high concentration 
(c > 1.5). Since high i nitial concentration should lead to 
a rapid collapse phase (IChernoff & Weinberglll990l) and no 
mechanism is known that could steepen the GMF over time, 
originally massive clusters with a shallow IMF and high c 
should still be visible in the "zone of avoidance" (upper right 
corner) of Figure[T] Admittedly, our sample is relatively small 
and we cannot exclude that clusters exist in this region of 
the parameter space. In particular, GCs of high concentra- 
tion (c > 2) and low total luminosity (My > -6) are potential 
candidates and their GMF should be investigated. 

Conversely, clusters with low concentration and steep IMF 
can have formed, NGC6809 being a good candidate. Objects 
like these should follow the general evolution towards lower 
a and possibly lower c. The exact balance between decrease 
in a and in c should depend on the initial conditions, and par- 
ticularly the initial mass (which together with the tidal radius, 
set by the Galactic potential, defines the relaxation time). 

Therefore, until GCs are found in the "zone of avoidance" 
of Figure[U it seems plausible that opposite evolutionary paths 
exist in the c-vs-a plane for clusters born with different cen- 
tral concentration and/or on different orbits, but that the IMF 
was the same or very similar for all of them. At low masses, 
the IMF of GCs must approach the steepest GMFs in our sam- 
ple, while at higher masses it cannot be much shallower than 
Salpeter without mass loss via stellar evolution causing the 
rapid disruption of the cluster (jChernoff & Weinberg! 1 19901: 
Takahashi & Portegies Zwart 2000). We presently have no di- 
rect measurements of the GC IMF above 0.8 M Q , but a value 
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of a ~ -2 is the preferred outcome of multi-mass Mic hie- 
King models dPDMt iDMPPt iDe Marchi & Pulonell2007h . In 
this sense, the tapered power-law dis t ributi on proposed by 
IDe Marchi. Paresce & Por tegies Zwart (2005) remains a vi- 
able hypothesis for the IMF of GCs. 

4. CONCLUSIONS 

We have discovered an empirical correlation between the 
central concentration and the GMF slope of GCs, whereby 
only loose clusters have a shallow GMF. A low value of 
the central concentration seems, therefore, a necessary con- 
dition for extensive mass loss leading to cluster disruption. 
Although it is possible that GCs formed with a certain de- 
gree of mass segregation and that some low-mass stars may 
have been lost due to tidal truncation before two-body re- 
laxation could act upon them, all of the depleted clusters 
that we have studied have a dynamical structure consis- 
tent with their being i n a cond i tion of energy equipartition 
dAndreuzzietal.1 120011: iDMPPt IDe Marchi & Pulonel 120071; 
iPulone. De Marchi & Parescell2007l) . Therefore, the observed 
trend between a and c can only be understood if either the 
depleted clusters have undergone collapse and have subse- 
quently rebounded, or if they are proceeding unnoticed to- 
wards core collapse or have already reached it without show- 
ing it. In either case, this means that the central concentration 
parameter c derived from the surface brightness profile is not 
a good tracer of a cluster's true central density or dynamical 
state. Our current estimate of the fraction of post-core col- 
lapse clusters may therefore need a complete revision as a 
large number of them may be lurking in the Milky Way. A 
reliable assessment of a cluster's dynamical state requires the 
study of the complete radial variation of its stellar MF. 
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